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Abstract 


In this paper, the general problem of a shell containing a 
through crack in one of the principal planes of curvature and 
under general skew- symmetric loading is considered. By employ- 
ing a Reissner type shell theory which takes into account the 
effect of transverse shear strains, all boundary conditions on 
the crack surfaces are satisfied separately. Consequently, 
unlike those obtained from the classical shell theory, the angu- 
lar distributions of the stress components around the crack tips 
are shown to be identical to the distributions obtained from the 
plane and anti-plane elasticity solutions. Extensive results are 
given for axially and circumferentially cracked cylindrical shells, 
spherical shells, and toroidal shells under uniform in-plane 
shearing, out of plane shearing, and torsion. Taking advantage 
of the fact that the problon is formulated for "specially" ortho- 
tropic materials, the effect of orthotropy on the results is 
also studied in some detail. 

1. Introduction 

In recent years, a great deal of effort has been devoted to 
the study of cracked plates and shells. The reason for this 
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steins from the fact that a great variety of advanced structures 
such as aerospace vehicles, pipelines, different components of 
nuclear reactors, etc. are designed from the viewpoint Of frac- 
ture mechanics. The early studies of cracked shells employed 
the so-called "classical" shell theory. This is an eighth order 
theory, and consequently can accommodate only four boundary con- 
ditions on each crack surface. To make the number of unknown 
functions arising from the solution of the differential equations 
compatible with the number of independent boundary conditions, 
the transverse shear and the twisting moment are combined as the 
"effective transverse shear" which in turn is used to satisfy the 
boundary condition regarding the transverse shear, (See, for 
example [1-7], and for review and references [8]). The short- 
comings of the classical shell theory are well-known from the 
plate and shell solutions. The angular distributions of the 
moments around the crack tips do not conform to those obtained 
from the plane elasticity solution and the transverse shear has 
a strong singularity of order of -3/2 which is not physically 
acceptable. It has been shown that in plates, if one uses a 
Reissner type plate theory [9-12] which adequately incorporates 
the effect of transverse shear strains, then the angular dis- 
tributions become identical to those obtained from the plane 
and anti-plane elasticity solutions [13-18]. Motivated by the 
success achieved in solving the plate problems, in recent years 
a Reissner type shell theory [19-20] has also been applied to 
cracked shells [21-26]. However, in all these studies, the 
loading is symmetric (i.e., the cracked shells are either 
under tension or bending).. In practice skew- symmetric 
loading is also very important. Considering also the fact 
that the existing solutions are limited, the need for a 
refined solution and additional reliable results does not 
need elaboration. In this paper the general problem of cracked. 
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thin and shallow shells under skew-symmetric loading will he 
formulated by using a Reissner type shell theory. After 
analyzing the asymptotic behavior of the stress resultants 
around the crack tips, extensive results for circumferentially 
and axially cracked cylinders, and spherical and toroidal shells 
under general skew-symmetric loading will be given. The effect 
of material orthotropy on the results will also be studied. 

2. Formulation of the Problem 

The problem under consideration is that of a shell of 
thickness h containing a crack of length 2a in one of the princi- 
pal planes of curvature and subjected to skew-synmetri c loading. 
The problem is formulated under the following assumptions: a) the 

shell is shallow and thin, b) the effect of transverse shear 
strains is included, and c) the material is "specially" ortho- 
tropic, i.e. the elastic constants are related as follows: 

^ C2.1) 

where 

E = /E-j E2 , V ~ 

Condition (.2.1) will enable us to factorize the differential 
equations. Referring to figure 1 for the geometry of the shell 
and notation, the equilibrium equations can be written as: 


ij,j ~ ° 

(2,2) 

i,i ^^*i ^-j^,j ° 

(.2.3) 

'ij,j-Vi=^ (i,j = l,2) 

(.2.4) 
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where the indicia! notation and the summation convention are 
used, and N-fj, M^j, V^- (i,j=l,2) denote the stress, 
moment and transverse shear resultants, respectively, q(X-|,X 2 ) 
is the load normal to the X-j-X 2 plane, and Z(Xi,X 2 ) is the equa- 
tion of the shell. Let U^*, W and 6^ (.i=l,2) he respectively the 
displacement components and the angles of rotation of the normal 
to the shell surface. Using the generalized Hook's law, in terms 
of the displacement derivatives, the strains are obtained as: 


=■ ^IjklNkl/h ■ ? W,^] (1.j = 1,2) 

C2.5) 


Defining a stress function F by 


Nlj=*1k=Jl f.kl (1.J-1.2) 


( 2 . 6 ) 


e||^ being the permuation symbol, the equilibrium eqs. (2.2) 
are satisfied identically. Eliminating U] and U 2 and substituting 
(2.6) into (2.5), we obtain: 


e- e. 6i e-i a* m tF* + h Z« . -e*! -■ = 0 

im jn kp Iq ijkr *mnpq ik jl *kl 


(2.7) 


For "specially" orthotropic materials satisfying condition (2.1), 
the strain-stress resultant relations may be expressed as: 


1 


^11 “ hE (N-n/c2-vN22)> g-j2 " Iie'^ ^12’ ^22 ~ hE (c2^22“'’’^ll ^ 


_ 1 


^11 " '^^‘^^^l,l'^^2,2^’ ^12 ^ ^^1,2'^^2,1^ 


( 2 . 8 ) 
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M22 - D(v6^^^+32,2/c ) » D - i2(l-v^) > 

(.2,9) 

chB ■ “,1 ■*■ ’ hB ^2 ^*2 ■*■ ® " 12(l+v} 

(2.10) 

where c = (Ei/E 2 )^. 


Using the relations (2.8)-(2.10) , the normalized 
defined in Appendix A, the curvatures defined by 

quantities 

3^Z _ 1 _ 1 32z _ 1 

^2 * " *^12 ’ 

(2.11) 

and the following new functions. 


36^ 33y 

n(x.y) " ly- - It > 

(2.12) 

33„ 33^ 

Kx.y) = k(— +^) - w , 

(2.13) 

equations (2.7), (2.3) and (2.4) become: 


V-+ - -JZ 3jz - 2 X,2^ 5),5y + Xj^ ^)w = 0 

(2.14) 

. X^d-KV^lCx,^ ^ - 2 x,2^ 3^ . X^a p , 


= xHI-i^v^) a 

(2.15) 

jC V^if; - i(; - W = 0 

(2.16) 

isi2^vaa-n = 0 

(2.17) 
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The problem is thus reduced to the solution pf the differential 
equations (2.14) - (2.17) under the boundary conditions as yet 
to be specified. 

For the crack problem the transverse load q=0, and noting 
that the crack is along a plane of principal curvature the 
shell parameter x -|2 = 0. Further it is assumed that through a 
proper superposition^ the problem is reduced to a perturbation 
problem where the only external loads are those acting on the 
crack surfaces; 

Defining 


V? = 


^1^ 


8 ^ 




3 ^ 


(2.18) 


eqs. (2.14) and (2.15) yield: . 

+ v2 (1 -kv2)((. = 0. (2.19) 


Now assuming (p in the following form, 

CO 

4>(Tr,y) = g(x,a)e"^^“ da 


( 2 . 20 ) 


and using the regularity condition at x=+», (2.19) gives: 


g(x,a) = 


f 4 m-x 

E R,-(a)e , X > 0 

1 J 


( 2 . 21 ) 


8 

E 

I 5 


m.x 


E R.-(a)e ^ , X < 0 

5 J 


Rj» (j=l».**»8) are unknown functions and mi,...,mg are the 
roots of the following characteristic equation: 
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= p + 

p** - K X 2 ** P^ + ( 2 KX-j^X 2 ^a^ ■ 2 icX 2 '*a^ + X 2 *^)p^ 

+ (2icXi^X2^a^ “ kX2^o^ ~ icX^'^a^ + 2X^ “ 2X^^X2^)a^p 

+ (X 22 - X^2)2 «*» = 0, 
such that, 

Re(mj) < 0 , nij +4 = -iHj j = 1,...,4 


w(x,y) = ^ I f(x,a)e"’^“ da 


00 

n(x,y) = ^ I h(x,a)e"’^“ da 

—00 

00 

'p(x,y) “ ^ I e(x,a)e da 


and substituting into (2.15)-(2.17) we obtain: 


f(x,a) = 


tf Ri(a)p ^2 m.x 




e ^ , X > 0 


, 8 Ri*(«)P?. 4 m-x 

X^ z T-2ST2 i 2^ e , X < 
c Xo ni-'^-X,'^a'^ 


5 A2 ...j -A^ 
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h(x,a) = 


A^(a) e , X > 0 
-r,x 

AgCa) e , X < 0 


(2.28) 


and 


e(x,a) = 


4 R,-(a)p,.2 m.x 

e , X > 0 




, 8 Ri(“)Pi-4 


(2.29) 
0 


where 


= - /■ 




(2.30) 


Considering the symmetry of loading and geometry with respect 
to the y-axis, the stress resultants satisfy the following sym- 
metry conditions: 

Nxx(x,y) = -N^j^(-x,y) , Nj^yCx.y) = N^^(-x,y) 

M^x(x.y) =? Mj^(-x.y) 

Vx(x.y) = V^(-x,y) (2.31) 


In solving the problem, it is therefore sufficient to consider 
the x>0 portion of the shell only. The boundary conditions of 
the problem can then be expressed as follows: 
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Nxx(0*y) =0 -« < y < 

00 


(2.32) 

Mxx(0»y) =0 -~ < y < 

00 


(2.33) 

lim N (x,y) = F,(y) 
x^+ ' 

1 

A 

A 

t ■ 

(2.34) 

v(0,y) =0 

|y| > J 



lim. M (x,y) = F 2 (y) 
x.^+ .xy ^ 

-i/c < y < i/c 

1 

> 

(2.35) 

3y(0,y) = 0 

\ 

A 



lim^ V (x,y) = F 3 (y) 
x->0+ ^ 

- t/q < y < 1 

1 

(2.36) 

w(0,y) = 0 

|y| > c 




where F-|(y), F 2 (y), F^Cy) define the loading on the crack sur- 
faces and are known functions. Using etjs. (2.6), (2.8)-(2.10) , 
(2.12), (2.13) and (2.20)-(2.30) the stress resultants may 
be expressed as: 

00 

if 4 in«x • 

Nxx(x»y) = “ ^ J J Rj(a)e ^ e"’^“ da (2.37) 

—00 

Nyy(x.y) = ^ I J dcs (2.38) 

-CO 

4 f ■ 4 ‘ m*x j.. 

Nj^(x,y) = ^ J a z mj Rj(a)e ^ e"’^“ da (2.39) 
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Application of boundary conditions (2.32)-(2.36) would lead to a 
system dual integral equations. However, if one defines a set 
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of new functions in terms of displacement and rotation deriva- 
tives, the problem can also be reduced to the solution of three 
simultaneous singular integral equations. 

Let, 


G^(y) = lim |^ - lim (^) y 


•5 

x-;b'+ 


o p 

GgCy) ~ lim. 


G3(y) = lim. 1^ 
^ x-^+ 


(2,45) 

(.2.46) 

(2,47) 


-lya 

e da 


Using the basic equations resulting from the formulation of the 
problem, the new unknown functions may be expressed as: 

^1^^^ \ (2.48) 

— 00 

, If" - 

' CO • 

+ K I a r-| A.j(a)e da (2.49) 

—00 

_• (■" 4 R. (a)p 4,, 

- " ■ 2?J “ f 

-00 ^ J • 

The homogeneous boundary conditions (2.32)-(2.33) and the 
inversion of (2.48)-(2.50) give: 


E R^(a) = 0 

1 J 


(2.51) 
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R,-{a)P,-^(ni4^-va2) 


' ■* •" “ ^l/^l 


= 0 


4 R,-(a)P,-^ 


(2.52) 

(2.53) 


X^ Jj Rj(a)p^^ 

T^ 1 (X2^m>-X^=^a^)(icPj-l) 

4 

1 = ^l(«) 

where 

qi(a) = I Gi(t)e’“* dt 

-./c 

q2(«) = j'^ G2(t)e^*“* dt 
-/c 

q3(a) = G3(t) e'“^ dt . 

-✓c 


2(ct) (2.54) 

(2.55) 

(2.56) 

(2.57) 

(2.58) 


By solving the system of equations (2.51)-(2.55) one can easily 
determine the unknown functions R-|(a),.. .,R^(a), A-j(a) in terms 
of G-|(y), G 2 (y), G 3 (y). The remaining mixed boundary conditions 
(2.34)-(2.36) give the following equations: 


lim^ 



4 m.x 

E m.R.(a)e ^ 

1 J J 



da = F^(y) 


✓c < y < i/c 


(2.59) 
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,00 




x-*0 


X2 


m.x 


4 R,-(a)p,-2m. , 


lya 


da 


Tim 'X a .. .0- v)2 


x->0+ 2ir hX^ 


I A,(c)(r,2+a2)e''^’' da = 


-/c < y < 


(2.60) 


Xr>2 it '"If 

x-^O"^ 2 tt , (Xg^nij^-Xi^a^jlKPj-l) 


e-y“ da 


- lim 
x->0 


i K “ A^(“)e ^ e da = f^Cy) 


-}/c < y < }/c 


(2.61) 


The integrands in (2.59)-(2.61) are bounded for every value of 
a except when |a[-»«> and y->t. Therefore one must extract the 
singular parts of the integrals when |a|-^. The functions R-j(a), 
..., R^Ca), A-jCa) may be written in the following form: 

Rj(a) = i[Qj(a)q-j(a) + Nj(a)q2(a) + Mj(a)q2(a)] 

(j=l,...,4) (2.62) 

A-|(a) = B^(a)q-|(a) + B 2 (a)q 2 (a) + B 2 (a)q 2 (a) (2.63) 

The functions Q. (a), N.(a), M.(a), (j=l,...,4) and B.(a), 

J J J * 

(i=l,2,3) are known from the solution of the system of equations 
(2.51)-(2.55). Now using the asymptotic expansion of m. , 
(j=l>--.»4) and r-| for large values of [a] 
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(2.64) 


nij(a) = “let 1(1 + 2 ^ ■ gjir + ...) J (j=l»...»4) 
ri(a) = “la|(l ” •••) > (2.65) 


after some lengthy analysis, the singular parts of the inte- 
grals which appear in (2.59)-(2.61 ) can be separated and we 
obtain: 


Vc G,(t) 3 |•/c 

-f-— dt + E k, .(y,t)6.(t)dt 

^ j-i U ' " 


= 2irF^(y) , -v^ < y < v'c 


(2.66) 


■I_y2 ^^(t) 3 

' ^ dt + E 

j = l 


j-VO U2V1 


= 2. f F2(y) 


fVE 

k2j(y.t)Gj(t)dt 

/c 

-t/c < y < Vc 


(2.67) 


y/c G-(t) 3 ri/c 

dt + E k3.(y,t)G.(t)dt 

i/E * U ' ' 


= 27rF3(y), -y/c < y < y/c . (2.68) 

The kernels k..(y,t) which appear in equations (2.66)-(2.68) 

I J 

are given in Appendix B. 

To complete the formulation of the problem, one must also 
impose the single-valuedness conditions for the displacements 
and for the rotation of the normal, which are: 
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9 


G2(t)dt-= 0 
G3(t)dt = 0 , 

-v^ 


rv^ 


P [G^(t) + (^) t G3(t)]dt = 0 


(2.69) 

(2.70) 

(2.71) 


The problem is thus reduced to the solution of the singular 
integral equations (2.66)-(2.68) under the single-valuedness 
conditions (2.69)-(2.71). 


3. Solution of the integral equations 

Once the system of integral equations (2.66)-(2.71) is 
solved, then all field quantities can easily be computed in 
terms of the functions G-|(y), G 2 (y)» G 3 (y). To do this, first 


the equations are normalized by defining: 

t = i/c^ T , < t < »^, -1 < T < 1 , (3.1) 
y = »^ n » <y</c,-l<n<l » (3.2) 
x = »^C,0<5,x<» (3.3) 
G.(»^t) = H.(t) , (i = 1,2,3) (3.4) 


Therefore, 

(3.5) 


The functions Hy(t), (i=l,2,3) are singular at t=+1 

h. (t) 

H.(t) = 1 R.(i= 1 . 2 , 3 ) 

’ (T+l)“(l-T)^ 
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where h^. ( t) are bounded and Holder-continuous in -1 < x < 1. 
Noting that the index of the problem k = C.a+e) = 1, and using 
the function-theoretic method described in [27] and [28], we 
have a = B ~ J • Then, the integral equations can be solved 
numerically by using a Gaussian quadrature type formula [29-31]. 

4. Asymptotic behavior around the crack tips. 

As stated earlier, one of the main objectives of this 
study is to show that by using a shell theory which adequately 
takes into account the effect of transverse shear strains, one 
can remove the discrepancy that exists between the classical 
shell theory and the elasticity solutions regarding the angular 
distribution of the stresses. Using the expressions given by 
(2.64) and (2.65), and the relation [28] 

I d. = exp[1(6 - 1^)] 

+ h(-l )exp[-i(B - ^ylj-)] + 0(-j-^)> (1 b|-^) (4.1) 

around n=l and c=0, the asymptotic expressions for the stress 
resultants are found to be: 

"xx ” f cos[B(l-n) - Bde (4.2) 

0 

, h^(l) -Bid _ 

N y = — ^ — cos [B(l-n) - x]dB (4.3) 
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”xy 


h,0) 


— sin[e(l-Ti) - J ]d6 

/3 


Ch^(l) 

2^ 


0 

/$ Sin[3(l-n) - J]d3 


M, 


XX 


a 

W 


(l-v2) f /3e'^l^l cos[3(l-n) 

J 


2'/^ 


I 

4 


M 


yy " hx^ 


hjO) fe-elsl 


= — 0-v2) I cosC6(l-n) - J]ds, 

0 

) p . r 

_ _ I " e-inro 

a ^ J ^ 


M = 1 h 2 

M„.. = T7r - — 


r®® - p ^ 

^ sin[6(1-n) - j]de 

J JK ^ 


1 h 


V.. = 


12 a 


*^tT 


2v^ 


00 

I ^ 


.-B|5| 


sin[3(l-n) - f]d3. 


0 

f3U| 


sin[3(l-n) - J]d3 , 


^ h.d) f -3UI 

v„ = -4rr cos[s(l-n) - Bde 

^ 'PH J /3 

0 


If we now define the new coordinates. 


? = r sine , n-1 = r cose 


and use the following expression [32] 


(4.4) 

■]d3 

(4.5) 

(4.6) 

(4.7) 

(4.8) 

(4.9) 

(4.10) 
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00 


^„-l ^-SZ. {S^,(rz)dz = . {=’"}(„ tan-’ 

cos ( 52 +y. 2 )lJ/'^ COS 

) 

F ’ 

(s > 0 , u > 0 ) , 

(4.11) 

(4.2)-(4.9) become: 


"xx = - C- J 1 * 4 2®] • 

(4,12) 

V *■’ ^ 2 4 sin 2 ] . 

(4.13) 

"xy = - 2 ^ f 4 1 ^ 4 2® ] • 

(4.14) 

M > h r 1 6 . 1 e-!n 50 

^xx 12a 2 ^ *- 4 2 4 2 ^ ’ 

(4.15) 

M z h ^2^^V r 7 0 1 50 n 

V ' 12 a 2 ^ *■ 4 2 4 sin 2 3 » 

(4.16) 

M- — h r 3 0 , 1 50 T 

^xy . 12 a 2 ^ [ 4 cos 2 + 4 cos ^ ] , 

(4.17) 

- a 

(4.18) 

- M’l n 

V = sin p . 

y 2 

(4.19) 


It should first be noted that the new coordinates defined by 
(4.10) are not polar coordinates for "specially" orthotropic 
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materials. For isotropic materials it is seen that the dis- 
tributions given by (4,12-4.19) are identical to those obtained 
from the plane stress and anti-plane elasticity solutions [18], 
[33]. 


5. Results and Discussion 

After solving the system of equations (2.66)-(2.71), the 
discrete values of the unknowns G-j(t), GgCt), GjCt) can be used 
to determine any desired field quantity in the shell. The 
mode n and mode IH stress intensity factors are defined as; 



Using expressions (4.12)-(4.19) and those given in Appendix A 
with 0 = 0 , the stress intensity factors can easily be expressed 
in terms of the end values of the unknown functions. In the 
examples given, three types of loading are considered. 

a) In-plane uniform shear loading: 

If denotes the magnitude of uniformly distributed 
shear stresses through the thickness, the external loads can 
be expressed as: • : 

r . - • 

Ni2(Q, X 2 ) = -Ni2 = » 


Xg) = 0 , 

V^(0, X 2 ) = 0, -a < X 2 < a [5.3) 


rl9- 



Using the dimensionless quantities defined in Appendix A, the 
functions F-j, F 2 , F^ that appear in (2.66)-(2.68) are found 
to be: 

F^(v^n) = - ^ , P2(.^ii) = 0 . n) = 0, 


-1 < n < 1 . (5.4) 

In this case the stress intensity factors will be normalized 
with . Thus, we have: 


_ 2 


k,(0) 


hid) E 


run 


m 


"m 


"tm 


_ k2(f)-k2(0) h2(l) ^ ^ 


cr i/iT 
m 


2a » 
m 


1^0 (0) 3R r— 

K«, = - 4 — hJl) 

m 


(5.5) 

(5.6) 

(5.7) 


b) Uniform Twisting Moment 

Let a^. denote the maximum of linearly distributed shear 
stresses through the thickness. Then, the external loads are: 

N^2^°’ ^ 2 ) = 0 , 

M-j2(0, X 2 ) = -M^2 " V '^t ’ 

V^(0, X 2 ) = 0 , -a < X 2 < a , (5.8) 

or 
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1 


F-jC*^ n) = 0, FgC/c i) = - -g » F^C*^ n) = 0 


-1 < n < 1 


(5.9) 


Normalizing the stress intensity factors with respect to » 
we obtain: 


kg(0) h,(1) g 




2 a. 


(5.10) 


"tt 


^St 


kg(h/2)-k3(0) _ h2(l)^j^ 




2 2a 


ko(0) O D 

-i— . -If v^hjd) 

o*/«r ^ “t ^ 


(5.11) 


(5.12) 


c) Uniform Transverse Shear Loading: 

If denotes the maximum of the parabolical ly distributed 
transverse shear stresses, the external loads can be expressed 
as: 

» t . 

! 

N^2(0> ^2^ " ° » 

Ml 2 Co* ^2) ~ 0 , 

v^(0, X2) = -v^ = - I h Og , -a < X2 < a (5,13) 


or 

FiCVc n) = 0, F2 (Vc n) = 0, F2(»'^ n) 


3 B ^ • 


-1 < n < 1 . 


(5.14) 


- 21 - 



In this case, normalizing the stress intensity factors with 
respect to , we obtain: 


_ - !lllV E 


ms 


05 *^ 


2 ’ 


'ts 


_ kg(h/2)-kg(0) hg(1) ^ ^ 


s 


2 03 2a ‘ 


(5.15) 


(5.16) 


1^^ 3 ^ 

” Oj/? 2 Oj 3' ' 


C5.17! 


The stress intensity factor ratios are calculated for different 
shell geometries. The following examples are considered; 

a) a cylindrical shell with an axial crack, 

b) a circumferentially cracked cylindrical shell, 

c) a spherical shell with a meridional crack, and 

d) a toroidal shell containing a crack at different locations. 

In all these examples, the material is assumed to be isotropic, 
with v=0.3. To give also some idea about the effect of ortho- 
tropy on the results, an example is considered for an axially 
cracked cylinder. The elastic properties of titanium used in 
these calculations are as follows: 

= 1.039 X 10^ Vm2 , 

E 2 = 1.434 X 10^^N/m2 , 

= 0.1966, 

V 2 = 0.2714, 
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5 


G^2 " ^-^75 X 10^° N/ni2 
Gave =4.956 x 10^° N/m2 , 

The results are given 1n tables 1-34. Some of the results are 
also displayed in figures 3 and 4. First, it must be pointed 
out that for toroidal shells, according to the crack location 
or X 2 ^ are taken as negative to describe negative curva- 
tures. Secondly, as x-j or X 2 approach zero, the results given 
in the tables approach those obtained for flat plates £183. 

From tables 1-34 and figures 3-4 the following trends can be 
observed: 

- If the crack surfaces are loaded with in-plane shearing 
forces, the coupling stress intensity factor ratios k^„, and 
*^sm small compared to As X 2 increases, i.e. as the 
curvature increases, k^ also increases. On the other hand, 

null 

a/h which expresses the thickness effect, does not affect the 
results significantly. 

- Under uniform twisting moment applied to the crack sur- 
faces, the stress intensity factor ratio k^^ becomes dominant, 
and the other two components k^^^^ and k^^ are small in compari- 
son. In this case, however, k^^ does not vary significantly with 
X 2 » and takes smaller values as a/h increases. 

- For uniform transverse shear loading, the stress intensity 
factor ratio k^^^ is very small, however, k^^ is no longer 
negligible. For this case also, the dominant stress intensity 
factor ratio is k^^. The stress intensity factor ratio k^^ 

does not vary significantly with curvature however, and decreases 
as the thickness of the shell increases. 
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- Figure 3 shows the comparison of for three shell 
geometries, assumes its smallest values for a cylinder 
with an axial crack, and assumes larger values for a circum- 
ferentially cracked cylindrical shell. The results for a 
spherical shell are even larger. 

- The results given in tables 28-31 for toroidal shells show 
the same trends for the stress intensity factors as the cylin- 
drical shells. It must be pointed out that, for this case, the 
component k^^ remains almost unaffected of the location of the 
crack. 

- The effect of orthotropy on the results is given in 

tables 33-34 and figure 4. As it meo' be observed, k^^^ is the 
component that is affected most by material orthotropy. If the 
axes of orthotropy are subjected to a 90® rotation, the results 
change significantly. For example, for E-|/E 2 < 1 are 
smaller than those obtained for isotropic materials, whereas 
for E 1 /E 2 ^ larger than the isotropic results. 
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Appendix A 


Dimensionless quantities used in the derivation: 


, X, X« X. 


u, 1 % U 

u- V-— W-- 


Px = , By = ^ 62 » •** ** 


a^ 


‘^xx = ^ * V = ‘'12/^ 


‘^xz " » V ' *^23/^ 


Nxx = N^^/chE , Nyy = cNgg/hE , 


xy 


”xx " M^i/ch^E , Myy = CM22/H2E , 
Vx = h B , V =Vc V2/h B 


= 12(l-v2) . V " ■•2(1 -v 2) ->h2R-2 


^12 " ‘•20-v^) = 12(l-v2) ^ ^ 


_ E 

» K 


(A.l) 

CA.2) 


(A. 3) 

(A.4) 
(A. 5) 
(A.6) 
(A. 7) 
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Appendix B 


Expressions of the kernels which appear in the integral 

4 

k-jiCy»t) = '1 [2ia Z nij.Qj(a)+l]sina(t-y)da 

O 

00 

" "I ^ mjNj(a)sina(t-y)da 
0 

00 ^ 

ki3(y,t) = -eJ a E mjMj(a)cosa(t-y)da 
0 

f" 1 . 4 Q4(a)P4^ni. 

kjity.t) = 2i j C-jr" J 


“ K (a^+ri^)B-j (a)]sina(t-y)da 


^22 


f” 1.4 N,(a)p.2m. 

(y.t) = j [21 ^ a r 


T) 


- i (o^+r^2)B2(a) - ■ ^ - ^~4^^] sina(t-y)da 

f”i_, 4 M.(a)p.2m. 

k23(y.t) = 2 J {rjz- « ^ (x2^m^.^-X^V)(icp.-l ) 

- (a2+r^2)B^(^^a)]cosa(t-y)da 


equations; 

(B.l) 

CB.2) 

CB.3) 

CB.4) 

(B,5) 

CB.6) 
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-i ►: a B^Ca)]cosa(t-y)da 


f 4 N.Ca)p^^ni. 

k32(y»t) = 2 J [i x2 K E (X zm Lx^^a^j^Kp.-l) 
O ^ J » J 

- 1 K a B2(a)]C0SaCt-y)da 

r 4 M.(a)p.3m. 

= 2 j [-X K S (xgi^mj^-xijVjcKPj.-l) 


+ K o BjCa) - 1 ]sina(t-y)da 


CB,71 


(B.8) 


CB.9) 
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Table 1. Stress intensity factor ratios k(nm for an iso- 
tropic cylinder containing an axial crack, v = 0 . 3 . 
(Uniform in- plane shearing) 


^1 

a / h=l 

a / h =2 

a / h =3 

- 


0.0 

1.000 

1.000 

1.000 

1.000 

0.25 

1.003 

1.003 

1.003 

1.003 

0.50 

1.012 

1.011 

1.011 

1.011 

0.75 

1.026 

1.024 

1.024 

1.024 

1.0 


1.040 

1.039 

1.039 

1.5 


1.078 

1.076 

1.075 




1.114 

1.113 

3.0 




1.182 


Table 2. Stress intensity factor ratios k^,^ for an iso- 
tropic cylinder containing an axial crack, v = 0 . 3 . 
(Uniform in-plane shearing) 



a / h=l 

a / h =2 

a / h =3 

a / h =4 

0.0 

0.000 

0.000 

0.000 

0.000 

0.25 

- 0.005 

- 0.002 

- 0.001 

- 0.001 

0.50 

- 0.012 


- 0.002 

0.000 

0.75 

- 0.018 

- 0.007 

- 0.001 

0.003 

1.0 


- 0.006 

0.003 

0.008 

1.5 


0.002 

0.017 

0.027 

2.0 



0.037 

0.051 

3.0 




0.103 
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Table 3 . Stress intensity factor ratios ksm for an isotropic 
cylinder containing an axial crack, v = 0 . 3 . (Uniform 
in-plane shearing) 



a / h=l 

a / h =2 

a / h =3 

a / h =4 

0.0 

0.000 

0.000 

0.000 

0.000 

0.25 

0.005 

0.004 

0.003 

0,003 

0.50 

0.020 

0.014 

0.012 

0.011 

0.75 

0.042 

0.030 

0.026 

0.024 

1.0 


0.051 

0.043 

0.039 

1.5 


0.101 

0.086 

0.078 

2.0 

j . 


0.137 

0.124 

t 

3.0 




0.222 . 


Table 4. Stress intensity factor ratios kmt for an isotropic 
cylinder containing an axial crack, v = 0 . 3 . (Uniform 
twisting moment) 



a / h=l 

a / h =2 

a / h =3 

a / h =4 


0.000 

0.000 

0.000 

0.000 


- 0.001 

- 0.001 

- 0.001 

0.000 




- 0.002 

- 0.001 



- 0.004 

- 0.003 




- 0.005 

- 0.004 

- 0.003 



- 0.007 

- 0.005 

- 0.004 

2.0 



- 0.006 

- 0.005 

3.0 




- 0.005 
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Table 5. Stress intensity factor ratios for an 

isotropic cylinder containing an axial crack, 
v=0.3. (Uniform twisting moment) 



a/h=l 

a/h=2 

a/h=3 

II 

0.0 

0.523 

0.354 

0,274 

0.228 

0.25 

0.523 

0.354 

0.274 

0.228 

0.50 

0.522 

0.353 

0.274 

0.227 

0.75 


0.353 

0.274- 

0;227 

1.0 


0.353 

0.274 

0.227 

1.5 


0.352 

0.273 

0.227 

2.0 



0.273 

. 0.227 

3.0 




0.226 


Table 6. Stress intensity factor ratios kst for an isotropic 
cylinder containing an axial crack, v=0.3. (Uni- 
form twisting moment) 


m 

a/h=l 

a/h=2 

a/h=3 

a/h=4 

0.0 



nm^m 

-0.091 

0.25 


-0.092 

-0.094 

nmm 

0.50 

-0.070 


-0.094 

-0.091 

0.75 

-0.069 

-0.091 

-0.094 

-0.091 

1.0 


-0.091 

-0.093 

-0.091 

1.5 


-0.089 

-0.092 

-0.090 

2.0 




-0.089 

3.0 




-0.087 
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Table 7. Stress intensity factor ratios k^,s for an iso- 
tropic cylinder containing an axial crack, v = 0 . 3 . 
(Uniform transverse shear loading) 



a / h=l 

a / h =2 

a / h =3 

a / h =4 

0.0 

0.000 

0.000 

0.000 

0.000 

0.25 

- 0.005 

- 0.005 

- 0.006 

- 0.007 

0.50 


- 0.021 

- 0.024 

- 0.028 

0.75 


- 0.045 

- 0.052 

- 0.061 

1.0 


- 0.076 

- 0.088 

- 0.103 



- 0. 149 

1 - 0.174 

- 0.204 




[ - 0.264 

- 0.310 

3.0 




- 0.495 


Table 8. Stress intensity factor ratios k^g iso- 

tropic cylinder containing an axial crack, v = 0 . 3 . 
(Uniform transverse shear loading) 



a / h=l 

a / h =2 

a / h =3 

a / h =4 

0.0 

0.47 

1.26 

2.12 

3.03 

0.25 

0.47 

1.26 

2.12 

3,03 

0.50 

0.46 

1.26 

2.12 

3.03 

0.75 

0.46 

1.25 

2.11 

3.02 

1.0 


1.24 

2.10 

3.00 

MM 


1.19 

2.04 

2.93 

Di 



1.95 , 

2.81 

3.0 




2.52 
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Table 9. Stress intensity factor ratios kgg for an isotropic 
cylinder containing an axial crack, v =0.3. (Uni- 
form transverse shear loading) 


m 

a / h=l 

a / h =2 

a / h =3 

a / h “4 

0.0 

1.68 

• 2.34 

2.98 

3.62 

0.25 

1.68 

2.34 

2.98 

3.62 

0.50 

1.67 

2.34 

2.98 

3.62 

0.75 

1.66 

2.33 

. 2.98 

3.61 

1.0 


2.32 

2.96 

3.60 

1.5 


2.28 

2.92 

3.54 

2.0 



2.84 

3,46 

3.0 




3.24 


Table 10. Stress intensity factor ratios kmj, for a cir- 
cumferentially cracked cylindrical shell; 


\2 

a / h=l 

a / h =2 

a / h =3 

a / h =4 

0.0 

1.000 

1.000 

1.000 

1.000 

0.25 

1.003 

1.003 

1.003 

1.003 

0.50 

1.013 

1.012 

1.012 

1.012 

0.75 

1.031 

1.027 

1.026 

1.026 

1.0 


1.047 

1.045 

1.045 

1.5 


1.104 

1.097 

1.095 

2.0 



1.165 

1.159 

3.0 




1.320 
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Table 11. Stress intensity factor ratios for, a cir- 
cumfereatially cracked cylindrical shell. 


X2 

a / h=l 

a / h =2 

a / h =3 

tl 

0.0 

0.000 

0,000 

0.000 

0.000 

0.25 

- 0.005 

- 0.002 

0.000 

0.001 

0.50 

- 0.014 

- 0.003 


0.006 

0.75 

- 0.023 

- 0.002 


0.017 

1.0 


0.004 

0.023 

0.035 

n 


0.030 

0.068 

0.092 




0.139 

0.178 

3.0 




0.437 


Table 12. Stress intensity factor ratios ks^ for a cir- 
cumferentially cracked cylindrical shell. 


^2 

a / h=l 

a / h =2 

a / h =3 

a / h =4 

0.0 

0.000 

• 

0.000 

p.ooo 

0.000 

0.25 

0.014 

0.010 

0.008 

0.007 

0.50 

0.056 

0.038 

0.032 

IIIIIIIIIIQIIgl^ 

0.75 

0.125 

0.084 

0.070 

0.062 

1.0 


0.148 

0.123 

0.109 

1.5 


0.328 

0.272 

0.242 

2.0 



0.479 

0.426 

3.0 




0.949 
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Table 13. Stress intensity factor ratios for a 

circumferentially cracked cylindrical shell. 


m 

a/h=l 

a/h=2 

a/h=3 

• 

II 

0.0 

O.OQQ 

0.000 

0.000 

.0,000 

0.25 


-0.001 

“0.001 

nmum 

0.50 

-0.005 


: ^0.002, 

-0.002 

0.75 

-0.009 

-0.005 . 

-0.004 

-0.003 



-0.007 

-0.005 

-0.004 

1.5 


-0.012 

mBtmm 


2.0 



-0.011 

-0.008 

3.0 




-0.012 


Table 14. Stress intensity factor ratios k^j. for a cir- 
cumferentially cracked cylindrical shell. 


X2 

a/h=l 

a/h=2 . 

a/h=3 

a/h==4 

0.0 

0.523 

0.354 

0.274 

0.228 

0.25 

0.523 

0.354 

0.274 

0.228 


0.522 

0.353 

0.274 

0.227 


0.521 

0.353 

0.274 

0,227 

1.0 


0.352 

0.273 

0.227 

1.5 


0.351 

0.273 

0.226 

2.0 



0.272 

0.226 

3.0 




0.225 
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Table 15. Stress intensity factor ratios kst for a cir- 
cumferentially cracked cylindrical shell. 


^2 

a/h=l 

' 

a/h=2 

a/h=3 

a/h=4 

0.0 

-0.070 

-0.092 

-0.094 

-0.091 

0.25 

-0.070 

-0.092 

-0.094 

-0.091 

0.50 

-0.070 

-0.091 

-0.094 

-0.091 

0.75 

-0.070 

-0.091 

. -0.094 

-0.091 

1.0 


-0.091 


-0.091 



-0.091 

-0.093 

. -0.091 

m 



-0.093 

-0.091 

3.0 




-0.091 


Table 16. Stress intensity factor ratios k^e for a cir- 
cumferentially cracked cylindrical shell. 


A2 . 

a/h=l 

a/h*2 

a/h=3 

a/h“4 

0.0 

0.000 

O.QOO 

0.000 

O.QOO 

0.25 

-0.016 

-0.015 



0.50 

■iBIiH 

-0.061 

-0.068 1 

-0.078 

0.75 

-0.133 

-0.133 

-0.149 

-0.172 

1.0 



-0.251 

-0.291 

1.5 


-0.387 

-0.462 

-0.546 

2.0 



-0.595 

-0.722 

3.0 



























































Table 17. Stress intensity factor ratios kts for a cir- 
cumferentially cracked cylindrical shell. 


X2 

a/h=l 

a/h=2 

a/h=3 

a/h=4 

0.0 

0.47 

1.26 

2.12 

3.03 

0.25 

0.46 

1.26 

2.12 

3.03 

0.50 

0.46 

1.25 

2.11 

3.01 

0.75 

0.43 

1.22 

27ot 

2.97 

1.0 


1.15 

1.99 

2.86 

im 


0.92. 

1.67 

2.46 




1.24 

1.89 

3.0 




0.93 


Table 18. Stress intensity factor ratios kgs for a cir- 
cumferentially cracked cylindrical shell. 


X2 

a/h=l 

a/h=2 

a/h=3 

a/h=4 

0.0 

1.68 

2.34 

2.98 

3.62 

0.25 

1.67 

2,34 

2.98 

3.62 

0.50 

1.65 

2.33. 

2.97 

3.60 

0.75 

1.58 

2.28 

2.92 

3.55 

1.0 


2. 17 

2.81 

3.44 

Esm 


1.77 

2.40 

3,00 

2.0 



1.84 

2.35 

3.0 




1.24 
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Table 19. Stress intensity factor ratios kn^ for a spheri- 
cal shell with a meridional crack. 


X 2 

a / h=l 

a / h =2 

a / h =3 

a / h =4 

0.0 

1.000 

1.000 

1.000 

1,000 

0.25 

1.006 

1.006 

1.006 

1.006 

0.50 

1.026 

1.024 

1.023 

1.023 

0.75 

1.058 

1.051 

1.049 

1.048 

1.0 


1.086 

1.082 

1.081 

1.5 


1.172 

1 . 162 . 

1.158 

2.0 



. 1.252 

1.244 

3.0 

i 



1.418 


Table 20. Stress intensity factor ratios ktm for a spheri- 
cal shell with a meridional cracL 


^2 

a / h=l 

a / h =2 

a / h -3 

a / h ?4 

0.0 

0.000 

0.000 

O.OQQ 

0.000 

0.25 

- 0.009 

- 0.004 

- 0.002 

- 0.001 

0.50 

- 0.021 

- 0.005 

0.003 

0.008 

0.75 

- 0.030 

- 0.001 

0.016 

0.025 

1.0 


0.011 

0.037 

0.052 

1.5 


056 

0.105 

0.136 

2.0 



0.205 

0.254 

3.0 




0.581 
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Table 21. Stress intensity factor ratios ksm for a spheri- 
cal shell with a meridional crack. 



a / h=l 

a / h =2 

a / h =3 

a / h =4 

0.0 

0.000 

0.000 

0.000 

0.000 

0.25 

0.014 

0.010 

0.009 

. 0.008 

0.50 

0.077 

0.053 

0.045 . 

. 0.040 

0.75 

0.169 

0.115 

0.096 

0.086 

1.0 


0.200 

0.166 

. 0.149 

1.5 


0.436 

0.361 . 

0.322 

2.0 



0.626 

0.558 

3.0 




1.209 


Table 22. Stress intensity factor ratios k^t for a spheri- 
cal shell with a meridional crack. 


m 

a / h=l 

^^9 

a / h =3 

a / li =4 

0.0 

0.000 

0.000 

0.000 

0.000 

0.25 

- 0.003 

- 0.002 

- 0.001 

- 0.001 

0.50 

- 0.008 

- 0.004 

- 0.003 

- 0.002 

0.75 

- 0.012 

- 0 . 007 . 

- 0.005 

- 0.004 

1.0 


- 0.010 

- 0.007 

- 0.005 

1.5 


- rO .014 



2.0 




- 0.009 

3.0 
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Table 23. Stress intensity factor ratios ktt Tor a spheri-. 
cal shell with a meridional crack. 


X2 

a/h=l 

a/h=2 

a/h=3 

a/h=4 

0.0 

0.523 

0.354 

0.274 . 

, 0.228 

0.25 

0.523 

0.354 

0.274 

0.227 

0.50 

0.521 

0.353 

0.274 

0.227 

0.75 

0.519 

0.352 

0.273 

0.227 

1.0 


0.351 

0.273 

0.227 

1.5 


0.350 

0.272 

.0.226 

2.0 



. 0.271 

. 0.226 

3.0, 




0.226 


Table 24. Stress intensity factor ratios kst for a spheri- 
cal shell with a meridional crack. 


\2 

a/h=l 

a/h=2 

a/h=3 

a/h=4 

■BH 


-0.092 

-0,094 

-0.091 

^J7Z5 

mmygiiim 


-0.094 

-0.091 

IIHiB 

-0.070 

-0,091 

-0.094 

-0.091 

^JT75 


-Tr:'09i ■ 

-0.093 

-0.091 

■mmn 


-0.091 

-0.093 " 

-0.091 

^175 


-0.090 

-0.093 

-O.09l 

“ZTD 



-0.092 






-0.090 
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Table 25. Stress intensity factor ratios k^is for a spheri 
cal shell with a meridional crack. 


X2 

a/h=l 

a/h=2 

a/h=3 

a/h.=4 

0.0 

0.000 

0.000 

0.000 

0,000 

0.25 

-0.023 

-0.022 

-0.025 


0.50 

-0.085 

-0.083 

-0.094 

-0.108 

0.75 



-0.199 

-0.230 



-0.282 

-0.325 

. -0.379 

1.5 


-0.465 

-0.563 

-0.669 

2.0 



-0.691 

-0.844 

3.0 




-0.843 


Table 26. Stress intensity factor ratios k^s for a spheri 
cal shell with a meridional crack. 


\2 

a/h=l 

a/h=2 

a/h=3 

a/h=4 

0.0 

0.47 

1.26 

2.12 

3,03 

0.25 

0.46 

1.26 

2.12 

3.03 

0.50 

0.45 

1.25 

2.11 

3.00 

0.75 

0.42 

1.20 

2.05 

2.94 

1.0 


1.12 

1.95 

2.81 

1.5 


0.86 

1.58 

2.35 

2.0 



1.15 


3.0 
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Table 27. Stress intensity factor ratios k^s for a spheri- 
cal shell with a meridional crack. 


\2 

a/h=l 

mim 

a/h=3 


0.0 

1.68 

2.34 

2.98 

3.62 

0.25 

1.67 

2.34 

2.98 

3.62 

0.50 

1.64 

2.32 

2.96 

3.60 

0.75 

1.55 

2.25 

2.90 


1.0 


2.12 . . 

2.77 

. 3.39 

1.5 


1.70 

2.32 

. 2.90 

2.0 



1.75. 

. 2.25 

3.0 






Table 28. The stress intensity factor ratios for the crack 
configuration shown in figure (2. a). Toroidal 
isotropic shell, with v=0.3. 



a/h=l 

Xi=0.575 

X2=0.257 

a/h=2 

Xi=1.150 

X2=0.514 

a/h=3 

Xi=1.725 

X2=0.771 

a/h=4 

Xi=2.300 

X2=1.028 

k 

mm 

1.019 

1.063 

1.116 

1.170 




0.047 

0.117 


0.040 

0.103 

0.180 


ft 

-0.006 

-0.007 

-0.007 

- -0.007 

^tt 

0.521 

0.352 

0.273 

0.227 

^st 

-0.069 

-0.090 

-0.092 

-0.089 

^ms 

-0.044 

-0.156 

-0.345 

-0.596 


0.46 

1.21 

1.93 

2.54 

^ss 

1.67 

2.28 

2.79 . 

3.18 . 






























































Table 29. The str«s intensity factor ratios for the crack 
configuration shown in figure (2.b). Toroidal 
isotropic shell, with v=0.3. 


1 

a/h=l 

Xl=0.257 

\2=0-575 

a/h=2 , 

Xl-0,514 

A2'=1,150 

a/h=3 

Al=0.771 

X2=1-725 

a/h=4 

Xl =1.028 
X2=2,300 


1.021 

1.073 

1,146 

1.231 


-0.019 

0.011 

0,107 

0,269 


0.079 

0.208 

0.384 

0.604 

l^nt ■ 

. -0.007 

. ^^0.009 

-0.010 

-0.010 


0.521 

0.352 

0.272 

0.225 

E39H 

-0.070 

-0.091 

-0,093 

-0,091 


-0.087 

-0,291 

-0.561 

-0,802 

mm 

0.45 

1.09 

1.47 

1.53 

^SS 

1.63 

2.06 

2,14 

1,95 


Table 30. The stress intensity factor ratios for the crack 
configuration shown in figure (2.c). Toroidal 
isotropic shell, with v=0.3. 



a/h=l 

a/h=2 

a/h=3 

a/h=4 


Xl=0.575 

Xi=1.150 

Xi =1.725 

Xi=2.300 


X2=0*332 

X2=0.664 

X2=0,996 

X2=1.328 


1.012 

1.039 

1.071 

1.105 

^tm 

0.008 

0.006 

0.003 

0.007 



-O.ooO 

0.011 , 

0.032 



07Do3 

Oo2 

0.001 

^tt 

0.522 

032 “ 

6.274 

0.228 

l^St 

-0.069 

-0.090 

-0.091 

-0.088 


-0.001 

-0.OO6 

-0.029 

-0.084 

*^ts 

0.4S 

1722 

1.95 

2.55 

^ss 

1.6^ 

^ 2723 

2.81 

3.17 
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Table 31. The stress intensity factor ratios for the crack 
configuration shown in figure (2.d). Toroidal 
isotropic shell, with v=0.3. 



Table 32. The stress intensity factor ratios for the crack 
configuration shown in figure (2.e). Toroidal 
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Table 33. 


Stress intensity factor ratios for an orthotropic 
(titanium) axially cracked cylindrical shell. 
(a/h=3, Ei/E2=0.725, Xq= [l2(l-v^)]l/‘ta/»^) 


X 

0 

non 


k 

sm 


^tt 


k 

ms 


k 

ss 


1.000- 

0.000 • 

0.000 

0,000 

0,274 

-0,091 

0.000 

2.07 

2.92 

0.25 

1.002 

-0.001 

.0,003 

rO.OOl 

0,274 

-0,091 

-0,005 

2.07 

2.92 


liOlO 

-0.003 

0.011 

-0,002 

0.273 

-0.091 

-0,020 ■ 

2.07 

2.92 

0.75 

1.021 

smm 

0.022 

-0,003 

0,273 

-0,091 

-0.043 

IHHElSii 

2.92 

1.0 

1.034 

0.001 

0,037 

-0.004 

0,273 

-0,090 

-0.074 

2.05 

2.91 

1.5 

1.067 

0.012 

0,074 

-0,005 

0,273 

-0.090 

-0.149 

2.00 

2.87 


1.103 

0.029 

0.118 

-0.006 

0.272 

-0.088 

-0.230 

1.93 

2.81 

























































Table 34. Stress intensity factor ratios for an orthotropic 
(titanium) axially cracked cylindrical shell . 
(a/h=3, Ei /E2=1.380, [12(1- v ^)]1/‘+ a/v^) 


X 

0 

k 

mm ' 


k 

sm 


^tt . 


k 

ms 

k 

ts 

k 

ss 

0.0 

1.000 

-0.000 

0.000 

0.000 

0.260 

-0.098 

-0.000 

2.11 

3.08 

0.25 

1.003 

-0.002 

0.004 

-0,001 

0,260 

-0,098 

-0.007 

2.11 

3.08 

0.50 

1.013 

-0.002 

0.015 

-0,002 

0,260 

-0,098 

-0.026 

2.11 

3.07 

0.75 

1.027 

-0.000 

0.031 

-0.003 

0,259 

-0,098 

-0.056 

2.10 

3.07 

1.0 

1.045 

0.006 

0.051 

-0.004 i 

0,259 

-0,097 

-0.095 

2.08 

3.05 

1.5 

1.085 

0.023 

0.101 

-0.006 

0.259 

-0,096 

-0.184 

2,01 

3.00 

2.0 

1.126 

0.047 

0,158 

-0,006 

0,259 

-0,095 

-0,274 

1.91 

2.91 





























Ri =3R 
Ro =5R 
Ri = R 
R2 =5R 
•i/Ri'O.I 
’’/Rg^O.Oa 










Fig. 3. Comparison of the stress intensity factor ratios 
Kmm for cylindrical and spherical shells 
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